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Abstract 

Subject to suitable boundary conditions being imposed, sharp inequalities are 
obtained on integrals over a region Q of certain special quadratic functions /(E) 
where E(x) derives from a potential U(x). With E = VU it is known that such 
sharp inequalities can be obtained when /(E) is a quasiconvex function and when 
U satisfies affine boundary conditions (i.e., for some matrix D, U = Dx on dU). 
Here we allow for other boundary conditions and for fields E that involve deriva- 
tives of a variety orders of U. We also treat integrals over of special quadratic 
functions g{J) where J(x) satisfies a differential constraint involving derivatives 
with, possibly, a variety of orders. The results generalize an example of Kang and 
the author in three spatial dimensions where J(x) is a 3 x 3 matrix valued field sat- 
isfying V- J = 0. We also present an algorithm for generating extremal quasiconvex 
quadratic functions. 

Keywords: Quasiconvexity, (5*-convexity, NuU-Lagrangians 

1 Introduction 



The divergence theorem 

[ V-U(x)rfx= [ V-ndS, (1.1) 

Jq JdQ 

of Lagrange, Gauss , Ostrogradsky, and Green, is one of the most important theorems in 
mathematics, particularly in applied mathematics. (Here n is the outward normal to the 
surface dQ of the region Q). In one dimension it reduces to the fundamental theorem of 
calculus and by setting U = wV one obtains the multidimensional version of integration 
by parts: 

/ V(x)-Vwrfx= [ wV-udS- [ w;V-V(x)rfx, (1.2) 
Jn Jan Jn 
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which with V = Vf yields Green's first identity, 



/ ■ dx = wVv ■ ndS — / wAv{x.) dx. 
Jn Jan Jn 

Note that the left hand side of (11.11) can be reexpressed as 



TrErfx, 



where E = VU is the gradient of the potential U. 

This immediately leads to the question: for what functions / can one express 



:i.3) 



:i.4) 



/(E(x))rfx 



:i.5) 



in terms of boundary data, when E derives from a potential? Such functions are known as 
null-Lagrangians. One example in two- dimensions of a null-Lagrangian is the determinant 
of E with E = VU and U being a two-component vector: one has 



/(E) = detE 



dUi dU2 dUi dU2 d 



dxi 8x2 9x2 dxi dxi 



dU2 
dxo 
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dxi 



:i.6) 



and we directly see that the quantity on the right is a divergence. With 



E 



El 
E2 



VU 

u 



:i.7) 



another null-Lagrangian is 

/(E) = (E2 ■ E2) Tr(Ei 



+ 2E2 ■ [(E1E2)] = (U ■ U)V ■ U + 2U ■ [(VU)U] 
= V-[U(U-U)], 



:i.^ 



which again is a divergence. It has been shown by Ball, Currie and Olver [T] that any 
null-Lagrangian can be expressed as a divergence, so that the evaluation of (11. 5p 
reduces to an application of the divergence theorem. Necessary and essentially sufficient 
algebraic conditions to determine whether a (quadratic or nonquadratic) function is a 
null Lagrangian have been given by Murat [21 [3l H] (see also [5]). In the important case 
where E = VU, /(E) is a null Lagrangian if and only if it is a linear combination of the 
subdeterminants (minors) of any order of E. For references see the paper of Ball, Currie, 
and Olver, who also show that when E(x) = V^U(x) there are no new null Lagrangians 
/(E) beyond those obtained by applying the result for = 1 to the ^-component 
potential V''-^U(x). 

Instead of seeking equalities it makes sense to search for sharp inequalities and try to 
find functions / and boundary data for which one can obtain sharp bounds on integrals 
of the form (11. 5p . If E = VU and affine boundary conditions U = Dx on dVl are 
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imposed, then one possible value of E in the interior of Q is of course E = D"^ (we take 
the convention that VUi, VU2, ■ ■ ■ , VUi are the columns, not rows, of VU, so that E is 
a. dx ^ matrix), and functions / that satisfy the inequality 

/ /(E(x))f/x> / /(D^)rfx=|r]|/(D^) (1.9) 

are called quasiconvex and for such functions the inequality is obviously sharp. (Here 
denotes the volume of Q). There is a large literature on examples of quasiconvex 
functions which include convex and polyconvex functions. So far only the quadratic 
quasiconvex functions have been completely characterized: see Tartar and Murat [6l [71 E] • 
In particular, the quadratic function /(VU) is quasiconvex if and only if /(H) is non- 
negative for all rank-one d x i matrices H. One of contributions of this paper is to show 
that for quadratic quasiconvex functions that are sharply quasiconvex, sharp bounds on 
the integral (11.51) can be obtained for a wide variety of boundary conditions, and not 
just affine ones. We emphasize (see p. 26 of |9] and [TO]) that it follows from an example 
of Terpstrafllj and has shown more simply by Serre|12[ [13] that there are quadratic 
quasiconvex functions that are not the sum of a convex quadratic function and a null- 
Lagrangian. (Terpstra's results imply that such quadratic forms exist only if £ > 3 and 
d > 3). This indicates that, in general, the sharp bounds we obtain cannot be obtained 
using nuU-Lagrangians. 

We allow for more general fields E(x), namely those with m components -E'.r(x), r = 
1,2, ... ,m, that derive from some real or complex potential U(x), with i components 
f/i(x), . . . , t/^(x) through the equations 

e 

Er{^)=Y,LrM^), (1.10) 

9=1 

for r = 1, 2, . . . , m, where Lrq is the differential operator 

h=l ai,...,a/i = l 1 ^ n 

of order t in a space of dimension d with real or complex valued constant coefficients Arq 
and ^"g^'"'' , some of which may be zero. 

We will find that sharp inequalities on the integral (11. 5p can be obtained for certain 
boundary conditions when / is a sharply Q*-convex function that is quadratic. A Q*- 
convex function is defined as a function that satisfies 



(/(E)) > /((E)), (1.12) 

for all periodic functions E = LU that derive from a potential U(x) that is the sum of 
a polynomial U°(x) and a periodic potential U^(x) . Here the angular brackets denote 
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volume averages over the unit cell of periodicity. (By all periodic functions we also mean 
for all primitive unit cells of periodicity, including parallelopiped shaped ones). In the 
degenerate pointed out to me by Marc Briane (private communication), there 

may be potentials that give rise to periodic fields, but are not expressible as the sum of 
a periodic part and a polynomial part, e.g. with L = d/dxi the potential U = Xi cos(x2) 
is not so expressible. However such potentials are not of interest to us. 
The function / is sharply Q*-convex if in addition one has the equality 

(/(E)) = /((E)), (1.13) 

for some non-constant periodic function E = LU that derives from a potential U(x) 
that is the sum of a polynomial U°(x) and a periodic potential U^(x) . We call the 
fields E and its potential U(x) (5*-special fields. When the elements Aj.q vanish and 
the coefficients ^"g^''^'' are zero for all h ^ t then Q*-convexity may be equivalent to 
quasiconvexity, but otherwise they are not equivalent: see p] for a precise definition of 
quasiconvexity. 

Note that there could be a sequence of functions Ei(x), E2(x), E3(x), . . . , E„(x) . . . 
each having the same average value E°, but not converging to E°, such that 

hm (/(E„)) = /(E°) (1.14) 

n— ^oo 

in which case we would call / marginally Q*-convex, but not sharply Q*-convex unless 
there existed a field E such that the equality f ll.lSp held. In this paper we are not 
interested in marginally (5*-convex functions which are not sharply Q*-convex. 

An associated problem, which we are also interested in, is to obtain sharp bounds on 
integrals of the form 

%(J(x))rfx (1.15) 



n 



for fields J(x) with m real or complex components Jr(x), r = 1,2, ... ,m, that satisfy 
the differential constraints 



J2LlrM^) = (1.16) 



for g = 1, 2, . . . , £, where 



in which the bar denotes complex conjugation. Observe that the operators in (11. lip and 
(I1.17P are formal adjoint s. 

Again we will find that sharp inequalities on the integral f ll.lSp can be obtained for 
certain boundary conditions when / is a sharply Q*-convex function that is quadratic. 
In this setting a Q*-convex function g is defined as a function that satisfies 

(^(J)) > ^((J)), (1.18) 
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for all periodic functions J satisfying L^J = 0. The function is sharply (5*-convex if in 
addition one has the equality 

{9(1)) = gm, (1-19) 

for some non-constant periodic function J(x) satisfying = 0. We will call J(x) a 
(5*-special field. 

Presumably, the assumption that / and g are quadratic is not essential, but if they 
are not it becomes a difficult task to find ones which are sharply Q*-convex. Also there 
should be some generalization of the theory presented here to allow for (5*-special fields 
E and J(x) which are not periodic, but we do not investigate this here. 

The theory developed here generalizes an example given in |14j, reviewed in the next 
section, that itself stemmed from developments in the calculus of variations, the theory 
of topology optimization, and the theory of composites: see the books dSl [HI \TE\ . 
A key component of this is the Fourier space methods developed by Tartar and Murat 
[HI [TJ E] in their theory of compensated compactness for determining the quasiconvexity 
of quadratic forms. 



2 An example 

To make the general analysis easier to follow we first review an example given in |14] . 
which was used to obtain sharp estimates of the volume occupied by an inclusion in a 
body from electrical impedance tomography measurements made at the surface of the 
body. 

Let us consider in dimension ci = 3 a 3 x 3 real valued matrix valued field J(x) 
satisfying V ■ J = 0. (Thus its three columns, not rows, are each divergence free). 
Alternatively J(x) could be regarded as a 9 component vector which corresponds to the 
case m = 9 and t = 1 in f ll.l6p and f ll.l7p . We want to show, that for certain fiuxes 
q = J-^n at the boundary dfl one can obtain sharp lower bounds on the integral f ll.lSp 
when 

^(J) = Tr(J^) + Tr(J^J) - [Tr(J)]l (2.1) 

This special function, introduced in section 25.7 of [17], is known to be quasiconvex 
(which in this case is equivalent to Q* convexity), in the sense that the inequality 

{giJ))-gi{3))>0 (2.2) 

holds for all periodic functions J(x), where the angular brackets denote volume averages 
over the unit cell of periodicity. Following the ideas of Tartar and Murat [HI El E] the 
reason can be seen in Fourier space where the inequality (by Parseval's theorem, since 
the function g is real and quadratic) takes the form 

Y,g{Re[3{k)]) + g{lm[3{k)]) > 0, (2.3) 
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which holds if g{ii) is non-negative for all rank two real matrices H, and in particular for 
the matrices Re[J(k)] and Im[J(k)] which are at most rank two since V • J = implies 
k ■ J(k) = 0. To prove that g{ii) is non-negative for all rank two real matrices, one first 
notes that g{ii) is rotationally invariant in the sense that 

^?(R^HR) = gin) (2.4) 

holds for all rotations R. This implies it suffices to check non-negativity with matrices 
H such that [0, 0, 1]^ ■ H = 0, i.e. of the form 



hu 


hi2 




h22 














H = I /l21 h22 h23 I , (2.5) 

in which case 

g{U) = {hn - h22f + (/^i2 + h2if + hl^ + hl^ (2.6) 

is clearly non-negative and zero only when hu = h22, ^12 = —^21, and /113 = /123 = 0. 
More generally, g{ii) will be zero if and only if the rank two matrix H has matrix elements 
of the form 

Hij = ao{kikj - 6ijk ■ k) + f3oeijmkm, (2.7) 

for some constants and Pq, and for some vector k which is a vector such that k-H = 0, 
where eijm is the completely antisymmetric Levi-Civita tensor taking the value +1 when 
ijm is an even permutation of 123, —1 when it is an odd permutation, and otherwise. 
This establishes (12. 2 p and shows one has the equality 

{gil)) - gm = 0, (2.8) 

for fields J(x) of the form 



where the are elements of a constant matrix J" and a(x) and /3(x) are arbitrary 
periodic functions (in this case J(x) has Fourier components of the form (12. 7p ). This can 
be established directly by noting that 

^7(1) = Tr(i^i), (2.10) 

where 

§ = 1 + 1^ - Tr(J)I = E° + 2VVa, E° = J° + (J°)^ - Tr(J°)I, (2.11) 

and so (12. 8p follows by integration by parts: 

{gm = (Tr(I^^)) = Tr{(J°)^[jO + (J^)^ - Tr(jO)I]} + 2(Tr[(J)^VVa]) 

= g{i')-2{iy-lfWa])=g{{l)). (2.12) 
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Now suppose J(x) defined within Q satisfies the boundary condition q = J-^^n at the 
boundary dfl where q£ has components 

= + - ^''^'^ + '''-d^ ) ' ('-'^^ 

for some functions a(x) and /3(x) defined in the neighborhood of 9f2, and for some 
constants J°£. Then one can extend a(x) and /3(x) beyond this boundary, so that they 
are periodic (with a unit cell of periodicity C containing VL) and define a field J(x) given 
by ( 12.91) . One can also extend the function J(x) so that it is periodic equals the function 
J(x) outside Vt and all its periodic copies. The boundary conditions fl2.13p ensure that 
this periodic function satisfies V ■ J(x) in a weak sense and therefore (12.21) implies 

! ^7(J(x)) rfx + / ^7(J(x)) rfx > |C|(7((J)), (2.14) 
Jc\vi Jn 

where |C| is the volume of C. Also ( 12. 8 p implies 

/ (7(I(x)) rfx + /" gUi^)) rfx = \C\g{a)). (2.15) 
Jc\n Jn 

Furthermore since V ■ J = weakly, integration by parts implies that 

[ Jkidyi= I V'^Jif(ix= I XkUiZie dS = [ i^rfx, (2.16) 

Jc Jc ~[ (^^i JdC ~l Jc 

and thus it follows that (J) = (J). So subtracting fl2.15p from fl2.14p gives 

/ (7(J(x))rfx> / (?(I(x))dx. (2.17) 
Jn Jn 

Using fl2.10p and integrating by parts, allows us to evaluate the right hand side: 

[ g{l{x))dx= [ Tr[(J)^V(x^Eo + 2V«)] c^x = [ n ■ l{x^Eo + 2Va) dS 
Jn Jn Jan 

= / q ■ (x^Eo + 2Va) (2.18) 
Jan 

In summary, for boundary fluxes of the form fl2.13p we obtain the inequality 

/ ^f(J(x))rfx > / q ■ (x^Eo + 2Va) rf^, (2.19) 
Jn Jan 

which is sharp, being satisfied as an equality when J(x) = J(x) within Q. In fact there 
is a huge range of fields for which one has equality since one is free to change a(x) and 
/3(x) in the interior of Q so long as the right side of fl2.13p remains unchanged at the 
boundary dQ. 
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3 Conditions for (3*-convexity 

Here we solve the problem of characterizing those quadratic functions f or g which are 
(5*-convex. This is a straightforward extension of the ideas developed by Tartar and 
Murat P, El [8] in their theory of compensated compactness for characterizing quadratic 
quasiconvex functions. Since / and g are quadratic we can let 

m m 

/(E(x)) = E ■ SE = J]5^E,(x)^,,E,(x), 

r=l s=l 
m m 

(7(J(x)) = J ■ TJ = J] 5^J,(x)T,,J,(x). (3.1) 

r=l s=l 

where the S^s or T^^ are the real or complex valued elements of some Hermitian m x m 
matrix S or T. 

We first consider periodic functions E(x) and J(x) and potentials U(x) that can be 
expressed in the form 

U(x) = U°(x) + Ui(x), (3.2) 
where U^(x) is periodic with zero average value and U°(x) is a polynomial with elements 

t d 

Uq = Bg + Bgl^ ^Xa^Xa2 ■ ■ ■ ^a^, (3-3) 

h=l ai,...,a|j=l 

where the coefficients Bg and B^j^'"""^ are chosen so that 

E'r=Y.Lr,U'g{^) (3.4) 
3=1 

is a constant for r = 1,2, ... ,m. We let denote the vector space spanned by all 
constant fields E'' with elements expressible in the form (13. 4p . 
By expanding E(x), U^(x) and J(x) in a Fourier series 

E(x) = (E) + ^ E(k)e^''-^ 
Ui(x) = ^Ui(k)e'''" 
J(x) = (J) + 5^J(k)e^^- (3.5) 

k^O 

where the sum is over all k in the reciprocal lattice space we see that the differential 
constraints fll.lOp or (11.161) imply that for all k 7^ in the reciprocal lattice 

£ m 

E,(k) = ^L,,(k)f^(k), ^Lt,(k)J,(k) = 0, (3.6) 

q=l r=l 
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where ^^^(k) and i^Jr(k) are the matrix elements 

t d 

Lrqi}^) = Arg + ^rqh ''^ai^a2 • • • ' 

h=l ai,...,a/i = l 

t d 

lI{\,) = a;+y, E (3.7) 

h=l ai,...,afi=l 

of the m X m matrix L(k) and its adjoint L'''(k), defined for all k G M^. In other words, 
the differential constraints imply that E(k) is in the range of L(k) and that J(k) is in 
the null-space J7k of L'''(k), where again these spaces are defined for all k G M'' and not 
just those k in the reciprocal lattice. We let ri(k) denote the projection onto £^k and 
r2(k) denote the projection onto jTk. Since these spaces are orthogonal complements it 
follows that 

ri(k) + r2(k) = 1, ri(k)r2(k) = r2(k)ri(k) = 0, vkGR"^ (3.8) 

where I is the m x m identity matrix. 

By substituting the Fourier expansion (13. 5p for E(x) in (11.121) we see that (11.121) holds 
if and only if the expression 

(E . SE) - (E) ■S{E) = J2 E(k) • SE(k) = /(E(k)) (3.9) 

is non-negative. A necessary and sufficient condition for this to hold for all possible 
primitive unit cells is that for all non-zero k G M'' 

/(H) > for all H G S^, (3.10) 

or equivalently that 

ri(k)sri(k) vkGR^k^^o, (3.11) 

which is an algebraic condition which can be checked numerically, and in some cases 
analytically. 

If S is real (and hence symmetric), then (13. 9p reduces to 
(E ■ SE) - (E) ■ S(E) = ERe[E(k)] ■ SRe[E(k)] + Elm[E(k)] ■ SIm[E(k)]. (3.12) 

k^O k^^O 

If in addition the coefficients Arq are real and the coefficients ^"g^ "'" are real when h is 

even and purely imaginary when h is odd then L(k) is real and H G £^k if and only if the 
real and imaginary parts of H lie in S^. Thus, in this case, to guarantee fll.l2p it suffices 
that fIXTUD holds for real H, for all k 7^ 0. 
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Similarly we look for Hermitian matrices T which are Q*-convex in the sense that 
fll.lSp holds and a necessary and sufficient condition for this is that 

^(H) > for all H G Jk, (3.13) 

or equivalently that 

r2(k)Tr2(k) Vk e k ^ 0. (3.14) 

If T is real and the coefficients Arq are real and the coefficients A'^^'f^'"''^ are real when h is 
even and purely imaginary when h is odd, then it suffices to check f l3.13p holds for real 
H. 



4 Sharply Q*-convex quadratic functions and their 
associated (3*-special fields 

The Q*-convexity conditions f l3.1ip and f l3.14p are clearly satisfied when S and T are 
positive definite, but in this case one only has equality in f lS.lOp and f l3.13p when H = 0. 
For any k 7^ let iSk denote the subspace of all m dimensional vectors H G for which 
one has equality in f l3.10p . and let denote the subspace of all £ dimensional vectors G 
such that L(k)G G S^. We call / sharply Q*-coiavex if, for some k 7^ 0, iSk contains at 
least one non-zero vector H. Note that if H G 5k then 

H-riSriH = 0, (4.1) 

and since from (13. lip FiSFi is positive semi-definite we deduce that H must in fact be 
a null-vector of this matrix, which implies that 

FiSH = 0, I.e. SH G Jk, (4.2) 

and hence that 

SL(k)GGJk, VGGWk- (4.3) 

Associated with a sharply (5*-convex / are periodic potentials U^(x) expressible in the 
form 

U^(x) = $^^(k)e*^■^ (4.4) 

where U^(k) G Wk for all k 7^ 0. The unit cell of periodicity has to be chosen so the k in 
the reciprocal lattice include some k such that contains at least one non-zero vector. 
Let us introduce the companion Q*-special fields 

E(x) = E° + ^E(k)e^''-'', 

J(x) = SE° + 5^ J(k)e*''-", (4.5) 
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in which G so that = LU (x) for some polynomial potential U (x) and where 

E(k) = L(k)n^(k), J(k) = SE(k). (4.6) 
With U(x) = U°(x) + U^(x) these Q*-special fields satisfy 

J(x) = SE(x), E = LU, L^J = 0, (4.7) 

and 

(E • SE) - (E) ■ S(E) = J2 ■ SE(k) = 0, (4.8) 

so that ( I1.13P is satisfied. 

Similarly for all k 7^ let 7k denote the subspace of all m dimensional vectors H G J7k 
for which one has equality in f l3.13p . We call g sharply Q*-convex if, for some k 7^ 0, 7k 
contains at least one non-zero vector H. Analogous to (14. 2p . (13.141) implies 

TsTH = 0, i.e. TH G S^- (4.9) 

Since £^k is the range of L(k), there exists a (possibly non-unique) £ dimensional vector 
G such that 

TH = L(k)G. (4.10) 
Now associated with a sharply Q*-convex g are the Q*-special fields 

I(x)=J° + ^2(kK''^ 

E(x) = +j2M'^y^"'' 

(4.11) 

where 

I(k) G Tk, |(k) = Tl(k) G £k, (4.12) 

and TJ° G so that TJ^ = LU°(x) for some polynomial potential U''(x). The unit 
cell of periodicity has to be chosen so the k in the reciprocal lattice include some k such 
that 7k contains at least one non-zero vector. We choose the periodic potential 

^i(x) = j^u^o^y^""^ (4-13) 

so that its Fourier coefficients satisfy 

|(k) = L(k)n'(k), (4.14) 
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for k ^ 0. With n(x) = n°(x) + ^^(x) these Q*-special fields satisfy 

^(x) = TJ(x), l = Ln, L'^J = 0, (4.15) 

and 

a ■ TJ) - U) ■ T{1) = ■ T2(k) = O, (4.16) 

so that f ll.iyp is satisfied. 

5 An algorithm for generating sharply (3*-convex quadratic 
functions and their associated (3*-special fields and 
for generating extremal quasiconvex functions 

To begin with let us suppose S can be expressed in the form 

S = V-vvVa, (5.1) 

where V is Hermitian and positive definite, and a is real and positive. The Q*-convexity 
of / is equivalent to the inequality 

H-VH>7(H-v)+7(lFV)-a|7|2 (5.2) 

holding for all complex 7 and for all H G S^, for all non-zero k G M'^, as can be established 
by taking the maximum of the right hand side over 7 thereby recovering (I3.10p . Rewriting 
this inequahty as 

a|7|^>7(H-v)+7(ir^)-H-VH, (5.3) 

and taking the maximum of the right hand side over H E we see that an equivalent 
condition is that 

a > V • r(k)v, Vk G k ^ 0, (5.4) 
where, assuming the i x i matrix L^(k)VL(k) is non-singular for all k 7^ 0, 

r(k) = L(k)(Lt(k)VL(k))-iD(k). (5.5) 

If V = I then r(k) is equal to ri(k), the projection onto the range of L(k). More 
generally r(k) is defined to be the matrix such that for all vectors v, 

r(k)vG£k and r(k)v-VvGJk, (5.6) 

implying 

r(k) = ri(k)[ri(k)vri(k)]-iri(k), (5.7) 

where the inverse is to be taken on the space £^k- In the context of quasiconvexity, the 
condition (15.41) was first derived in [TH], stimulated by a result of Kohn and Lipton [2U] . 
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The simpler derivation presented here was suggested by an anonymous referee of a later 
paper (see also section 24.9 in [T7]). 
Clearly ( 15 ■4p is just satisfied if 

a = sup V • r(k)v, (5.8) 
k e M"' 


and in this case / is marginally Q*-coiwex. However if this supremum is attained for 
some k 7^ then / is sharply Q*-convex. 
As an example, consider fields of the form 

E^(7). (5,9) 

where U{x.) is a scalar potential. Suppose V = I and 

t 

1 

where t is a real (i-component unit vector. Then r(k) is the projection 

1 



(5.10) 



r(k) = r,(k) = -^(^j(k- 1), (5.11) 



in which /c^ = k ■ k and so 

v.r(k)v=(i^^ (5,12) 

takes its maximum value a = 2 when k = t. (It is clear that the maximum occurs when 
t and k are parallel and then it is just a matter of simple algebra to show fl5.12p is less 
than or equal to 2.) This example shows that the supremum can be attained, and in 
some cases may be attained only at one non-zero value of k. 

Returning to the general case, let /C denote the set of those k 7^ which attain the 
supremum in (15. 8p . Then for k G /C, a = v ■ r(k)v and the set consists of those G of 
the form 

G = a(D(k)VL(k))-iD(k)v, (5.13) 
where a is an arbitrary complex constant. Thus 

H = L(k)G = aL(k)(L"f(k)VL(k))-iD(k)v = aT{k)v (5.14) 

satisfies 

SH = VH - av = -aA(k)V-V, (5.15) 

where 

A(k) = V - Vr(k)V (5.16) 
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has the property that 

D(k)A(k) = D(k)V - D(k)VL(k)(D(k)VL(k))-i&(k)V = 0, (5.17) 

implying SH G J]^. So the potential U^(x) and the (5*-special fields E(x) and J(x) take 
the form 

Ui(x) = ^a(k)[(Lt(k)VL(k))-iLt(k)]ve-'i'-, 
E(x) = E° + ^ a(k)r(k)ve'''-^, 

J(x) = SE(x) = SE° - J] a(k) A(k) V- V e**^'", (5.18) 

where E° e £° and 

a(k) = if k ^ /C. (5.19) 

The reciprocal lattice needs to be chosen so that it includes some k G /C to ensure the 
a(k) in fl5.18p are not all zero. 

When V = I then A(k) is equal to the projection r2(k) onto the null-space of 
L^(k). More generally A(k) is defined to be the matrix such that for all vectors w, 

A(k)w G Jk and A(k)w - V^^w G S]^, (5.20) 

implying 

A(k) = r2(k)[r2(k)v-ir2(k)]-ir2(k), (5.21) 

where the inverse is to be taken on the space Jk. 

We now use the fact that for an appropriate choice of the scalar polynomial ^(k), 

P(k) = g(k)(Lt(k)VL(k))-^Lt(k)v (5.22) 

is also polynomial in k: to see this it suffices to take 

g(k) = det[Lt(k)VL(k)], (5.23) 

although it is better to take the lowest degree polynomial which works. Then setting 

a(k)=g(k)c(k), (5.24) 

where the c'(k) are the Fourier components of some scalar periodic function c(x), we see 
that 

U^(x) = ^{d/dxi, d/dx2, . . . , (9/(9xd)c(x), (5.25) 

which also gives the Q*-special fields E = E° + LU^ and J = SE. The function c(x) has 
to be chosen so that (:(k) is zero if k ^ /C. (We do not have to worry about this latter 
constraint if it happens that K, consists of all non-zero vectors in M°'). 
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We assumed V was positive definite but in fact it is easy to clieck tliat tlie derivation 
goes tfirougfi if the quadratic form associated with V is Q*-convex. It could also be 
marginally Q*-convex but in this case one needs to choose v so that 

Lt(k)v lies in the range of Lt(k)VL(k) Vk ^ 0, (5.26) 

since otherwise the supremum in (15. 8 p is surely infinite. Thus by induction, successively 
setting 

V = Vo - ^VivJ/oi, v = Vj, a = aj, (5.27) 

1=1 

for j = 1, 2, . . . , n (where Vq has an associated Q*-convex quadratic form, v = has to 
be chosen so that f l5.26p is satisfied, and a = aj is given by (15. 8p ) one finally reaches a 
point such that the quadratic function / associated with 

n 

S = Vo- J]v,vJM, (5.28) 

i=l 

is an extremal quadratic Q*-convex function in the sense that one cannot substract any 
non-zero positive semi-definite tensor from S and retain its Q*-convexity. In this case, 
with V = S the supremum in (15.81) is infinite for all non-zero v G C". Incidentally, this 
solves the problem of generating extremal quadratic quasiconvex functions (or equiva- 
lently "extremal translations"): see [19] and chapters 24 and 25 in [17j for other insights 
into this problem. 

In the end one obtains Q*-special potentials of the form 

n 

Ui(x) = J2 Pjid/dxu d/dx2, . . . , d/dxd)c,{^). (5.29) 
i=i 

Here the n scalar periodic potentials Cj(x) have Fourier components satisfying 

Cjik) = if k ^ /Cj, (5.30) 

and, when one has made the substitutions (I5.27p . the Pj are polynomials given by (I5.22p 
and the )Cj consist of those k which attain the supremum in (15.81) . 
Similarly, supposing T can be expressed in the form 

T = - wwV/3, (5.31) 

where V is positive definite, we see that g is Q*-convex if and only if 

/3 > w ■ A(k)w, Vk G ky^ 0. (5.32) 

If we choose 

/3= sup w-A(k)w (5.33) 
k G M'^ 
k^^O 
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then g will be marginally Q*-convex and sharply Q*-convex if this supremum is attained 
for some k 7^ 0. Let Ai denote the set of those k 7^ which attain the supremum. Then 
for k e Al, ^ = w • A(k)w and the set consists of those H of the form 

H - 6A(k)w, (5.34) 

where b is an arbitrary complex constant. In this case we have 

TH = 6V-^A(k)w -bw = -6r(k)Vw = L\k)G, (5.35) 

where 

G = -6(L+(k)VL(k))-^L^(k)Vw. (5.36) 
So the potential (x) and the Q* -special fields E(x) and J(x) take the form 

^i(x) = -5^6(k)[(L^(k)VL(k))-^Lt(k)]Vwe*-", 
§(x) = TJ(x) = TJ° - ^6(k)r(k)Vwe*'^-'', 

I(x) = J° + XI ^(k) A(k)w e*'^-^, (5.37) 

k^O 

where e £° and 

6(k) = if k ^ A^. (5.38) 
We next look for the polynomial q{k) of lowest degree such that 

P(k) = -g(k)(Lt(k)vL(k))-iLt(k)Vw (5.39) 

is also polynomial in k. Set 

b{k) = g(k)c(k), (5.40) 

where the cr(k) are the Fourier components of some scalar periodic function c(x) such 
that now c(k) is zero if k ^ . Then we see that 

n'(x) = P{d/dxi, d/dx2, d/dxd)c{-x), (5.41) 

which also gives the g*-special fields E = T J° + LU^ and J = T'^E (in the case where 
T is non-singular). 

The derivation still goes through if the quadratic form associated with is Q*- 
convex or marginally (5*-convex. In the latter case one needs to choose w so that 

r2(k)w lies in the range of r2(k)V-^r2(k) Vk 7^ 0, (5.42) 
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since otherwise, from f l5.2ip . the supremum in f l5.33p is surely infinite. Thus by induction, 
successively setting 



i-i 

V-i=Wo-5^w,wJ/A, w = w„ /3 = /3„ (5.43) 

i=l 

for j = 1,2, ... ,n (where Wq has an associated Q*-convex quadratic form, w = wj has 
to be chosen so that fl5.26p is satisfied, and /3 = (3j is given by f l5.33p ) one finally reaches 
a point such that the quadratic function g associated with 

n 

T = Wo-^w,wl/a, (5.44) 

i=l 

is an extremal quadratic Q*-convex function in the sense that one cannot substract any 
non-zero positive semi-definite tensor from T and retain its Q*-convexity. In this case, 
with = T the supremum in f l5.33p is infinite for all non-zero w G C". 

6 Sharp inequalities on the integrals over Q 

Since the operators L and are formal adjoints the quantities 

fo= f /(E(x)) rfx = / E(x) • SE(x) = / E(x) ■ J(x), 

J(x)-E(x) (6.1) 



n 



90= [ g{W))d^= I I(x)-SJ( 
Jn Jn 



n 



can be computed in terms of boundary terms using integration by parts. Now we show 
that 

/ /(E(x)) rfx > /o, / (?(J(x)) rfx > go, (6.2) 
Jn Jn 

for all fields E(x) deriving from a potential U(x) that matches the appropriate boundary 
data of the potential U(x) or for all fields J(x) that match the appropriate boundary 
data of J(x), provided certain further ancillary conditions hold. First note that these 
inequalities are clearly sharp, being attained when E(x) = E(x) or when J(x) = J(x). 

To establish the first inequality in ( 16. 2 p we first find a parallelopiped C that contains 
Q and that is formed from an integer number of the primitive unit cells of E(x). (If Q 
lies inside a primitive unit cell of E(x), then we can take C as this primitive cell, but 
otherwise we need to join a set of these primitive cells together to obtain a parallelopiped 
that covers fl). We extend E(x) outside fl so that it is periodic with C as a unit cell. 
In this cell, but outside Q, it E(x) equals E(x). Since the potential U(x) matches 
the boundary data of the potential U(x) at the boundary dQ, the appropriate jump 
conditions are satisfied across this interface. We extend the potential U(x) outside C so 
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that U^(x) = U(x) — U (x) is C-periodic: if xq is at the center of another unit cell then 
for X e C we set 

U(x + xo) = U(x) + U°(x + xo) - U°(x), (6.3) 

which in particular gives the potential inside the copies of Q outside C. Defined in this 
way the relation E = LU holds in a weak sense, and so the inequality f ll.l2p is satisfied, 
which we rewrite as 

/ /(E)rfx+ / /(E)rfx> |C|/((E)), (6.4) 
Jn Jc\n 

where |C| is the volume of C. Also (14. 8 p holds, which we rewrite as 

fo+ [ /(E)rfx=|C|/((E)), (6.5) 
Jc\n 

So subtracting these equations gives 

/ /(E)rfx>/o + |C|[/((E))-/((E))]. (6.6) 
Jn 

This inequality in general requires us to know (E) and (E). Given a m- dimensional 
constant vector J" we have 

[ J° ■ (E - E) dx = [ J° ■ L(U - U) rfx = [ J° ■ A(U - U) rfx, (6.7) 
Jc Jc Jc 

where A is the matrix with elements A^q. (In establishing the last equality in (16. 7p we 
have used integration by parts and the fact that U and U are equal in the vicinity of the 
boundary of C). Since this holds for all J° we deduce that 

(E) = (E) + ^^A(U-U)dx. (6.8) 

So a sufficient condition for /((E)) to be equal to /((E)) is that 

SA = 0, (6.9) 

so that the range of A is in the null space of S (which if S is non-singular requires that 
A = 0). When this ancillary condition holds then clearly (16. 6p implies the first inequality 
in (I62D 

To establish the second inequality in (16. 2p we first find a parallelopiped C that contains 
Q and that is formed from an integer number of the primitive unit cells of J(x). Then 
we extend J(x) outside Q so that it is C-periodic, and within the unit cell C equals J(x) 
outside fl. Since the field J(x) matches the appropriate boundary data of the field J(x) 
at the boundary dQ, the appropriate jump conditions are satisfied across this interface. 
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Defined in this way the relation L^J = holds in a weak sense, and so the inequality 
fll.lSp is satisfied, which we rewrite as 



/ ^(J) dx + / gil) > \C\g{{3)). (6.10) 
JQ Jc\n 

Also fl4.16p holds, which we rewrite as 

go+ [ gil)d^=\C\gi{J)), (6.11) 
Jc\n 

and subtracting these equations gives 

/ g{3)d^>go + \C\[g{{3))-g{{3))]. (6.12) 

This inequality in general requires us to know (J) and (J). Given a m dimensional 
constant vector E° G £^ so that E° = LU°(x) for some polynomial potential U''(x) we 
have 

[ (J - J) • E° rfx = [ (J - J) ■ LU° dx = 0, (6.13) 
Jc Jc 

where the last equality follows from integration by parts, using the fact that J = J in a 

vicinity of the boundary of C. Since this holds for all E° G we deduce that a sufficient 

condition for g{(E)) to be equal to g{(E)) is that the range of T be a subset of £^ and if 

this ancillary condition holds then clearly f l6.12p implies the second inequality in (16. 2p . 



7 A Generalization 

Following the ideas of [21] (see also sections 24.1 and 30.4 of [T7]) the arguments present 
here extend directly to integrals of the form 

/i(E(x), J(x))rfx, (7.1) 

where h is quadratic and the m component field J(x) still satisfies = while the 
n-component field E(x) derives from a p component potential U(x) according to the 
relations 

^ p ^ ^ 

E,(x) = ^L,,f/,(x), (7.2) 

9=1 

for r = 1, 2, . . . , n, where L^q is the differential operator 

h=l a^,...,an=l OXa^ UXa; OXa^^ 
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of order f in a space of dimension d with real or complex valued constant coefficients Arq 
and<--\ 

The function /?.(E, J) can be expressed in the form 

ft(E.J)=(f).M(f), (7.4) 

where the matrix M is Hermitian. It is chosen to be sharply (5*-convex so that for 
periodic fields E(x) and J(x) satisfying the differential constraints one has 

(ME,J))>M(E),(J)), (7.5) 

with equality for some non-constant (5*-special fields E = E and J = J. 
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